
function [tt , xx] = Bcr4bp_Orbit3D(coord , tVec, x0, aux)
%
% BCR4BP画轨道，3D
%
% 时间：2024年06月03日
% 作者：张晨
% 单位：CSU
% 邮箱：chenzhang@csu.ac.cn
%%%%%%%%%%%%%%%%%%%%%%%

options = odeset('Reltol', 1e-8, 'AbsTol', 1e-8) ;

if strcmp(coord , 'emrot')

    [tt , xx] = ode113(@Bcr4bp_EqmEMRot3D , [tVec(1) , tVec(2)] , x0 , options , aux);

elseif strcmp(coord , 'sb1rot')

    [tt , xx] = ode113(@Bcr4bp_EqmSB1Rot3D , [tVec(1) , tVec(2)] , x0 , options , aux);

else

    fprintf('动力学错误！')
    tt = [];
    xx = [];
    return;

end

end
